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Polarization dynamics of femtosecond light pulses propa- 
gating in air is studied by computer simulation. A rich variety 
of dynamics is found that depends on the initial polarization 
state and power of the pulse. Effects of polarization on the 
plasma and supercontinuum generation are also discussed. 



I. INTRODUCTION 

There has been a great interest in long distance fem- 
tosecond pulse propagation in air in recent years moti- 
vated in part by potential applications in remote sens- 
ing and laser-induced lightning. The first experimental 
observations p]|| of highly localized, high-intensity fila- 
ments propagating over distances that exceed their cor- 
responding Rayleigh lengths by orders of magnitude mo- 
tivated the efforts to understand the phenomenon |j] |) . 
Due to the violent formation process, and the extreme 
time and spatial scales of the filaments, the computer 
simulations and analytic approaches turned out to be vi- 
tal tools to grasp the underlying physics |9-16|. Several 
models have been proposed. The first explanation sug- 
gested a stationary waveguide formed by competing ef- 
fects of nonlinear self-focusing and linear defocusing by 
underdenseplasma generated in the most intense part of 
the pulse M. An alternative model employed a notion 
of the "moving focus" to explain how the geometrical 
focus of a beam is transformed into a long filament ||. 
More recently, a dynamic spatial replenishment model 
emerged from numerical investigation by Mlejnek et al. 
P,^0|! The essential feature of the femtosecond propa- 
gation of infrared (IR) pulses is its dynamics that makes 
it possible that the localized filaments propagate over 
long distances and do not suffer significant energy loses. 
The basic mechanism involves a dynamic balance be- 
tween the nonlinear self-focusing and defocusing by free 
electrons generated via multi-photon absorption by the 
high-intensity filaments. When the self-focusing leading 
edge of the pulse starts to generate plasma, it thus cre- 
ates a defocusing "lens" for the trailing portion of the 



pulse. This has two effects; First, it limits the losses due 
to absorption in the plasma and, second, it prevents the 
major portion of the pulse from experiencing self-focusing 
collapse. After the plasma generating leading portion of 
the pulse exhausts its energy, the strength of the defo- 
cusing lens decreases, and the self-focusing starts over 
again, this time in the "next temporal slice" of the pulse. 
The whole process can repeat several times, depending 
on the total energy of the pulse. The qualitative features 
of the dynamical spatial replenishment remain valid also 
in transversely wide pulses, that break-up into multiple 
filaments p2]. 

Up to the present, most of the work on the exper- 
imental side and all numerical studies have been con- 
cerned with the case of linearly polarized input pulses. 
Recently, Petit et al. (Tt]] studied the effects of the po- 
larization on the propagation of femtosecond IR pulses. 
They have measured luminescence from the plasma gen- 
erated in the filaments to show that the polarization of 
the pulse plays an important role in the plasma genera- 
tion. Due to the highly dynamic nature of the filament 
formation and propagation, it is natural to expect a rich 
polarization dynamics in femtosecond pulses. In this ini- 
tial study, we restrict ourselves to femtosecond pulses 
with modest peak powers that retain their initial radial 
transverse symmetry and are just sufficient to produce 
several refocusing events within a single pulse. The ini- 
tial polarization of the pulse is varied and the polarization 
state is recorded along the propagation path. Our results 
indicate a tight correlation between the evolution of the 
pulse waveform along the propagation distance and cer- 
tain global polarization parameters. Thus, the measure- 
ments of the polarization state could provide yet another 
tool to extract information on dynamics of pulses as well 
as an opportunity to correlate experiment and theory. 

The remainder of the paper is organized as follows. 
Our basic model is described in Sec. II, and details of 
our simulated experiment are given in Sec. III. Plasma 
production, polarization dynamics, and the associated 
supercontinuum generation are then discussed in Sees. 
IV- VI. Finally, our summary and conclusions are given 
in Sec. VII. 
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II. MODEL EQUATIONS 

Since the typical transverse dimension of the self- 
focusing collapsing filaments are a few tens to hundred 
microns in the IR wavelength region, it is a good ap- 
proximation to consider the optical field as transverse. 
We describe the complex optical amplitude in terms of 
two circularly polarized components, £^. The choice of 
the circularly polarized basis is motivated by the fact 
that in this basis the nonlinear birefringence is "diago- 
nal" and preserves the power in each component, which 
in turn makes its implementation easier. Our model is 
a straightforward extension of the scalar version we used 
in our previous studies p,^2[. It takes into account the 
effects of diffraction, normal group velocity dispersion 
(GVD), multi-photon ionization (MPI) and avalanche 
plasma generation, defocusing of light by the generated 
free electrons, and instantaneous and delayed cubic non- 
linearity: 
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Here uj is the optical frequency, |£| 2 = |£ + | 2 + |£~| 2 the 
combined intensity of the two circularly polarized com- 
ponents, k = to/c, k" = d 2 k/dui 2 , p is the electron den- 
sity, a is the cross-section for inverse bremsstrahlung, r is 
the electron collision time, f}( K > is the AT-photon absorp- 
tion coefficient, and the nonlinear change in refractive- 
index for a linearly polarized continuous wave (cw) field is 
n 2 |£| 2 . The corresponding critical power for self- focusing 
collapse for linearly polarized fields is P cr — X 2 /2im2, or 
P cr =1.7 GW for our parameters. In contrast, the crit- 
ical power for self-focusing for circularly polarized fields 
P^, = 1.5P cr is 1.5 times that for linearly polarized fields 
p8| . The normalized response function (characterized by 
the resonance frequency f2 and the decay T) 



R(t) = 0(t) n 2 e" rt/2 
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A = v/f! 2 - T 2 /4 , (2) 



accounts for delayed nonlinear effects, and / is the frac- 
tion of the cw nonlinear optical response which has its 
origin in the delayed component, and we denoted the 
Heaviside step function by 9(t). In the present model, 
we chose the relative weight of the "self" and "cross" 
non-linear birefringent terms the same as for the instan- 
taneous Kerr effect and nonlinear Raman effect, namely, 
the cross effect has a weight twice that of the self effect 
as is appropriate for an isotropic medium |18|. 

The optical field amplitude equations are completed by 
a simple equation describing the evolution of the plasma 
density: 
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where the first and second terms represent the avalanche 
and multi-photon effects, while the last one models the 
plasma recombination. Note, that at the time-scales rele- 
vant for the present study, the only practically important 
contribution comes from the the multi-photon term. For 
the pulse powers we use, it is sufficient to include oxygen 
alone as a source for MPI, since its corresponding multi- 
photon order is lower than that for nitrogen. We use the 
Keldysh theory formula to calculate the MPI rate . 

Explanations of symbols that appear in our model 
equations are listed in Table 1 together with the values 
used in our simulation. 



III. NUMERICAL SIMULATION 

To study the role of the initial polarization state of 
the pulse on its subsequent propagation, we performed 
a series of simulations. In all runs, the initial pulse was 
chosen to be a Gaussian wave-packet (in time and space) 
characterized through its central wavelength A = 775 
nm, pulse waist w = 0.7 mm and temporal duration 
tfwhm = 200 fs. We fixed the initial peak intensity as 
1.0 x 10 16 W m -2 which is a relatively modest value: for a 
linearly polarized input this corresponds to an input peak 
power of P = 7.7 GW= 4.5P cr , whereas for a pure circu- 
lar polarization P = 7.7 GW= 3.0F C ^. At these powers, 
there are typically two to three refocusing events in the 
pulse propagation, and the numerics can be reliably con- 
trolled. At higher powers under a perfect axial symmetry 
of the pulse, it may be necessary to extend the model be- 
yond the nonlinear Schrodinger equation (NLSE) , and to 
include correction terms that allow to handle pulses with 
a very broad spectrum. To check our numerical proce- 
dures and model implementation, we performed most of 
the runs in two or three different resolutions. The data 
we present were obtained with the time-domain resolu- 
tion of 0.24 fs. This resolution is sufficient to capture 
spectra several hundredth of nanometers wide. We pro- 
vide more details of our testing procedures concerning 
the spectral resolution and algorithm validity in Sec. VI 
devoted to supercontinuum generation. 

Below, we present our results for a series of runs that 
differ in the initial polarization state of the pulse. We 
change the polarization from the linear, through ellip- 
tic to circular, to see how it affects the dynamics of 
the filaments. In previous work on nonlinear propaga- 
tion in fibers, the Stokes parameter formalism has been 
employed to classify the polarization dynamics for plane 
wave fields Q . Here we employ space and time averaged 
Stokes parameters [sq, s\, S2, S3) as a description of the 
polarization dynamics of the propagating pulses. In par- 
ticular, our Stokes parameters are calculated numerically 
according to the prescription 
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.S = (J>/2 + -^o) 

Si = (J>/4 - J ?r _ 7r /4)/so 

S3 = (^cr - ^cl)/s (4) 

where is the total energy detected after passing the 
pulse through a polarizer of state (3 

/ + oo r-K 
\ \£ p \ 2 {r,t)rdrdt (5) 
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with Bp the field resolved along the polarizer direction, 
and R the radius of a detection aperture chosen to select 
the the most intense region of the pulse around the fila- 
ment. For a linear polarizer (3 is the angle of the polarizer, 
whereas for a circular polarizer (3 — CR,CL correspond- 
ing to right and left circular polarization settings. Here 
we chose R = 0.1 mm, and perform the measurement of 
the Stokes parameters "in the near field". We note that 
,s is the total energy of the pulse detected over the aper- 
ture, and the remaining Stokes parameters are calculated 
as differences between the detected energy with different 
polarizer settings normalized to the total energy. Thus 
our prescription has a direct experimental interpretation 
and should therefore be of utility. The above prescrip- 
tion for determining the Stokes parameters also reduces 
to the usual definitions in the limit of long pulses of broad 
transverse extent. 

In the numerical simulations to be presented we fix the 
initial Stokes parameter S2 equal to zero, and vary Si and 
S3 between zero and one to vary the initial polarization 
from linear through elliptic to pure circular polarization. 
Besides the polarization state, we also recorded the data 
pertaining to the plasma generation inside the pulse, and 
generation of the supercontinuum light. We start our 
discussion with plasma generation. 

IV. PLASMA PRODUCTION 

As the femtosecond pulse undergoes multiple self- 
focusings, the amount of plasma generated by its high- 
intensity portions reflects the spatio-temporal shape of 
the pulse. The total number of generated electrons as 
well as the maximal plasma densities exhibit peaks along 
the propagation distance. These peaks coincide with the 
locations of self-focusing collapses within the pulse, each 
peak being produced by a different temporal portion of 
the pulse. Figure [j] shows the plasma generation for three 
different initial polarizations of the pulse. The trend that 
one can see is quite in line with what is expected based 
on the functional form of the nonlinear birefringence. 
Namely, as we change the initial polarization from linear 
through elliptic to circular, the amount of the generated 
plasma decreases. Also, the onset of filament formation 
is delayed for the circularly polarized pulse because the 
critical power for self-focusing is higher for circularly po- 
larized pulse as noted earlier. In other words, keeping the 



input peak power the same for different polarizations, 
we effectively decrease the self-focusing power of circu- 
larly polarized pulses. This is also the reason why the 
number of refocusing events can be higher in a close-to- 
linear or linear polarization than in a circularly polarized 
pulse. While the overall plasma production depends on 
the polarization state, the typical dimensions of the fil- 
aments are not very sensitive to it. That can be seen 
from the Figure |l| which shows the longitudinal extent 
of the plasma columns. The transverse dimensions of 
the plasma channels can be estimated from the ratio of 
the two curves shown in the figure as the square root of 
the ratio between the total number of electrons and the 
maximal plasma density. This characteristic dimension 
of the plasma channel is shown in Fig. || for three differ- 
ent polarizations. Though there are small variation be- 
tween different initial polarizations, the thickness of the 
plasma channels is always roughly 60 microns in the most 
dense parts. The plasma channel generated by the circu- 
larly polarized pulse seems to be more "homogeneous", 
exhibiting less thickness variation along the propagation 
distance. 



V. POLARIZATION DYNAMICS 

While the dynamics of the plasma generation and its 
dependence on the polarization described in the previ- 
ous Section is straightforwardly linked to the structure 
of the equations governing the optical field evolution, the 
polarization dynamics seems to be more difficult to in- 
terpret. Figures ||,||,|] show the Stokes parameters and 
the polarization degree as functions of the propagation 
distance for the three different initial polarizations we 
discussed in the previous Section. An interesting feature 
is the difference between the "stability" of initial linear 
and circular polarization. Figure [j] shows that a small 
perturbation to the linear polarization in the initial pulse 
leads to an increasing deviation of the polarization state 
from the initial one. In this sense, the linear polarization 
appears to be unstable, as the polarization measured af- 
ter the filament formation can significantly differ from 
the initial one. Naturally, the rate of divergence for two 
close but not identical initial conditions decreases with 
the decreasing input power. On the other hand, in the 
case of almost circular polarization shown in Fig. the 
pulse polarization state doesn't change that dramatically. 
Though there is a small decrease of the polarization de- 
gree, one can say that final polarization stays close to the 
initial one even after two refocusings of the pulse. Thus, 
the circular polarization seems to be more stable against 
small perturbations than the linear polarization. Figure^ 
shows an interesting case of an initially elliptic polariza- 
tion. Note that the Stokes parameter S3, which measures 
the degree of circular polarization, only exhibits small 
variations, while the other two parameters decrease sig- 
nificantly after their initial increase in the first collapse. 
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That means that the light focused in the second collapse 
is predominantly circularly polarized. This observation is 
confirmed by examining the spatio-temporal polarization 
pattern within the pulse. This is an observation that may 
not be expected based on the previous results concerning 
cw self-focusing of polarized pulses pi|Jlq| . Namely, in 
a situation close to a continuous wave regime, one can 
argue that the weaker circular component experiences 
a stronger focusing "lens" because of the factor two in 
the birefringence cross-term, and that eventually leads to 
equal intensities of both circular components and, there- 
fore, linear polarization of the central filament. However, 
the important point here is that the femtosecond light 
filaments under consideration are extremely dynamic ob- 
jects. The resulting polarization distributions strongly 
depend on the spatial and temporal location within the 
pulse, and any interpretation based on steady-state-like 
considerations becomes invalid. Namely, there is a delay 
between when the light encounters the focusing "lens" 
and when it actually reaches the focus. This delay inter- 
feres with the temporal profile of the pulse, which typi- 
cally exhibits multiple peaks that may be just a few fem- 
tosecond long. As a consequence, the above simple ar- 
gument is not sufficient to capture all essential features 
of the phenomenon. The tendency of the predominantly 
circular polarization of the most intense portion of a fila- 
ment was also observed in our simulations that were not 
restricted to axial symmetry j2^|. A wide beam with a 
random perturbation breaks up into multiple filaments 
that exhibit polarization properties similar to those we 
discuss here. However, one has to keep in mind, that 
in both cases, axisymmetric as well as fully spatially 
resolved, our simulation modeled pulses with relatively 
small energy fluence when compared to some current ex- 
periments. Therefore, it would be extremely interesting 
to see what happens to the polarization of the central 
filament in a pulse that has enough energy for many self- 
focusing events and also has the transverse profile clean 
enough to preserve its axial symmetry. 

We conclude this Section with yet another presenta- 
tion of the polarization dynamics data we have shown 
above. Namely, we want to show that the polarization 
changes closely reflect the self-focusing events within the 
pulse and, consequently, the locations where most of the 
plasma is generated. Figure ^| shows the root-mean- 
square rate of the change of the Stokes vector along the 
propagation distance 




for the case of elliptic initial polarization of the pulse. 
The curve shown corresponds to the data depicted in 
Fig. [5] and in Fig. [l]b). Note, that the maxima of the 
rate of the polarization change closely follow those in 
the plasma production curve. We thus see that the mul- 
tiple self-focusing events in the single pulse leave their 
signature on the polarization. This could provide an- 



other way, besides the indirect plasma density observa- 
tions |p|,p3|-^5|, to visualize the dynamics of the spatial 
replenishment. 

VI. SUPERCONTINUUM GENERATION 

After contrasting the behavior of pulses polarized close 
to linear and circular from the points of view of plasma 
generation and of their polarization dynamics, we want 
to discuss the effects of polarization on the supercontin- 
uum generation. However, before presenting our results, 
we feel a note concerning some technical questions is in 
order. The explosive spectral broadening in the super- 
continuum generation in femtosecond pulses is a rather 
subtle phenomenon from the point of view of numerical 
simulation. Clearly, one needs a sufficient resolution in 
the time (spectral) domain to capture the broad spec- 
trum, but the resolution may not be the only issue here. 
One has to check how broad is the spectral region within 
which the model and its numerical implementation de- 
scribes the wave propagation correctly. It is expected, 
that at extreme powers, correction terms beyond the ba- 
sic NLSE (see e.g. Ref. |p6j ) need to be included in the 
field equation. To ensure that we work in the regime 
where the correction terms may be neglected, and to as- 
sess the spectral band over which our numerics works 
well, we performed some comparative simulations. The 
choice of the reference frequency (wavelength) around 
which the NLSE is built is in principle arbitrary, though 
it is obviously most appropriate to choose it close to the 
central frequency of the modeled pulse. This means, that 
simulations that only differ in the choice of the reference 
frequency should give the same results. We have com- 
pared simulation with the reference frequency shift of 150 
nanometers off the central wavelength of the pulse, and 
obtained a very close match of the spectra in the region 
from 500 to 1200 nm and over four decades in spectral 
intensity. Thus, in this interval we can trust the spectra 
extracted from our simulations. We would like to point 
out that this is in fact a rather strong test for the overall 
numerical implementation of the solver. But most impor- 
tant, the close agreement between the spectra shows that 
the correction terms beyond NLSE do not play an impor- 
tant role in our regime of modest input powers; This is 
because one can interpret them as corrections that par- 
tially restore the invariance of the original wave equation 
with respect to the choice of the (physically meaningless) 
reference frequency - that is exactly what our test shows 
is not needed for the conditions in our simulations. 

Figure [?] shows comparison of the spectral broadening 
of two pulses that differ only in their initial polarization. 
The way we extracted the spectra from the pulse wave- 
forms corresponds to measurement in a near field with 
the same aperture we used for polarization characteriza- 
tion. The figure shows spectra "measured" after the last 
sclf-focusing collapse, after which the pulse will eventu- 
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ally diffract and there will not be more supercontinuum 
generation. One can see that the pulse which was initially 
polarized close to linear exhibits a much stronger super- 
continuum generation. However, note that we compare 
pulses with the same peak power, and what we see here is 
an effect similar to the plasma generation. Supercontin- 
uum generation strongly depends on the available power 
and the natural measure of that power is in units of crit- 
ical power for self-focusing. From that point of view the 
circularly polarized pulse is weaker and that is the main 
reason that it exhibits less spectral broadening. However, 
this situation represents a reasonable experimental setup 
in which only the polarization is changed. 

To get a feeling regarding the role of the group ve- 
locity dispersion in the supercontinuum generation, wc 
performed most of the simulation runs also with a higher 
group velocity dispersion parameter. It turns out that in- 
creasing GVD by an order of magnitude leads to a strong 
suppression of the continuum production. We speculate 
that it may be one of the reasons that, at least in some 
experiments, there is only little spectral broadening in 
the ultraviolet (UV) femtosecond pulses |^7| because the 
GVD value of air is significantly higher in the UV region. 

The findings from our numerical simulations should 
be accessible to experimental testing. However, extreme 
caution should be exercised when trying to compare ex- 
perimental and simulational spectra. The spectra we 
present are taken from "one shot". They exhibit mod- 
ulations typical for supercontinuum generation in gases 
P,p8| . While the characteristic "frequency" of the mod- 
ulation is rather reproducible, the exact spectral shape 
is not. Thus, even small fluctuation in the parameters of 
the pulses will result in a suppression of fine features in 
the multiple-shot experimental spectra. 

VII. CONCLUSIONS 

We have performed a computer simulation study of 
the effects of the initial pulse polarization on its prop- 
agation and filamentation dynamics. In agreement with 
the experiment [jlTj , we have found that the filamentation 
onset is reached earlier for a linearly polarized pulse than 
in a circularly polarized pulse of the same peak power. 
However, in some cases, the experiment indicates that 
the circularly polarized pulses create higher plasma den- 
sities in comparison with linearly polarized pulses. Our 
simulations suggest the opposite, but one has to keep 
in mind that the experimental measurements and our 
simulation pertain to rather different conditions, includ- 
ing much higher power and focusing in the experiment. 
In our simulations, we also see more self-focusing events 
with linearly polarized pulses than with circularly polar- 
ized ones. 

We have observed a rich spatio-temporal polarization 
dynamics. Naturally, the limited range of the full param- 
eter space explored in our simulations prevents us from 



drawing general conclusions, but we believe some tenden- 
cies are already discernible. First, the initially circular 
polarization seems to be stable in a sense that a small 
perturbation of the polarization state doesn't grow sig- 
nificantly. On the other hand, when a small polarization 
perturbation is applied to the linearly polarized pulse, 
it grows and the polarization degree of the central fil- 
ament decreases significantly. The growth rate of the 
deviation is expected to increase for higher powers. An 
interesting case is the one of an initially elliptic polariza- 
tion. We have observed that the center of the filament is 
almost purely circularly polarized after subsequent self- 
focusing collapses within the pulse. Apparently, we have 
here a rather different situation than in the self-focusing 
in nanosecond pulses, which tend to create linearly polar- 
ized filaments independently of initial polarization state 
|pl| , ^8[ . The usual argument, that the the weaker circu- 
lar component experiences stronger self-focusing which fi- 
nally leads to balancing the power of the circularly polar- 
ized components and thus to linear polarization, can't be 
applied in the femtosecond pulses. The evolution of the 
polarization state along the propagation distance is ex- 
tremely dynamic and hardly possible to describe in sim- 
ple "static" terms. The multiple refocusing within a pulse 
and the defocusing effect of the generated plasma play a 
major role. The important feature of the whole process 
is that apart from the relatively small energy losses due 
to plasma generation and radiation, the energy in each 
circular component remains conserved. Thus, the main 
mechanism that results in changing polarization pattern 
within the pulse is the spatio-temporal energy redistribu- 
tion within each circular component. As a consequence, 
the polarization state of the whole pulse is very compli- 
cated and, therefore, any projection onto global quanti- 
ties like Stokes parameters has to be interpreted in re- 
lation to details of the measurement (aperture, near vs. 
far field, collecting angle, . . . ). 

We have also looked at spectra "measured" after the 
last self-focusing collapse for different initial polariza- 
tions. In accordance with our observation about the 
plasma generation, we see much stronger supercontin- 
uum generation in linearly polarized than in circularly 
polarized pulses. 

Finally, we have seen that rate of change of the polar- 
ization state in the center of the filament is closely cor- 
related with self-focusing and plasma generation. Thus, 
the polarization offers, in principle, an alternative way to 
investigate the dynamics of the spatial replenishment in 
femtosecond pulses. 

In this work, we have concentrated on investigating 
"global" quantities to characterize the femtosecond pulse 
propagation that should be experimentally accessible, at 
least in principle. Naturally, the question is how much 
can be done practically. To measure the Stokes parame- 
ters evolution along the propagation distance, for exam- 
ple, a very good reproducibility of the initial pulse would 
be required. However, the "final output" polarization 
state and spectra, which should be easier to measure, 
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also carry a lot of signatures about the inner dynamics 
of the femtosecond pulse propagation. 
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FIG. 1. Plasma generation for almost linear (a), elliptic 
(b) and a close-to-circular (c) initial pulse polarization. The 
total peak power is kept the same in all cases. Since the crit- 
ical power for self-focusing is higher for circular polarization, 
the circularly polarized pulse experiences weaker self-focusing 
which in turn results in less overall plasma generation. 
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FIG. 2. Characteristic transverse dimension of the gen- 
erated plasma column for three different initial polarizations. 
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FIG. 3. Stokes polarization parameters of the on-axis part 
of the pulse as functions of the propagation distance for an 
almost linear initial polarization. The initial deviation from 
the perfect linear polarization increases, and the polarization 
degree decreases. 
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FIG. 4. An initially almost circularly polarized pulse (ap- 
proximately) preserves its polarization state. There is only a 
slight decrease of the polarization degree. 



10 



CO 

£ 

oS 




to 

Cl 

CO 
CD 

O 

CO 



c 



-0.5 



-1 



S2 



S1 
S2 
S3 

polarization degree 







1 



2 3 
propagation distance [m] 

FIG. 5. Stokes polarization parameters as functions of the 
propagation distance for an elliptic initial polarization of the 
pulse. The central part of the filament evolves into a predomi- 
nantly circular polarization state after the second self-focusing 
collapse event. Note that the rate of change of the polarization 
state correlates with the loci of maximal plasma production 
(see Fig. |). 



11 



8.0 




propagation distance [m] 

FIG. 6. Root-mean-square of the rate of change of the 
Stokes polarization vector, ds/dz (in units of m _1 ) as a func- 
tion of the propagation distance for elliptic initial polariza- 
tion. The rate maxima are correlated with the locations of 
strongest plasma generation and focusing. 
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FIG. 7. Spectral broadening of two femtosecond pulses 
with different initial polarizations. The linearly polarized 
pulse produces significantly more supercontinuum light than 
an equally intense circularly polarized pulse. 
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TABLE I. Model parameters and numerical values used in 
our simulations. 



Quant. 


Value&Unit 


Note 


k 


k = 2tt/\ 


reference wavevector 


Ao 


775 x 10 m 


wavelength 


k 


o i ~i n — 29 2 / 

2.1 x 10 s /m 


group velocity dispersion 




5.6 x 10 cm /W 


nonlinear index 


f 


0.5 




r 


26 THz 


R(t) ~ 6(t)x 


A 


16 THz 


e 1 t/2 sin(At) 




7 


MPT order 


E g 


fa 11 eV 


ionization energy 


13k 


6.5 x 10- 104 m n W- 6 


MPI rate 


T 


3.5 x 10 -13 s 


electron collision time 


a 


5 x 10~ 24 m 2 


cross-section for 






inverse bremsstrahlung 


a 


5 x 10~ 13 m 3 /s 


recombination rate 
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